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I. INTRODUCTION 

In order to understand the infrared properties of Quantum Chromodynamics [QCD], the heavy quark sector is 
a useful area of study. Among the heavy quark correlations, the most basic quantity is the confinement potential 
between a quark-antiquark pair. At large separations, Wilson loops in lattice calculations exhibit an area law which 
corresponds to a linearly rising potential whose coefficient, the Wilsonian string tension, can be explicitly related to a 
hadronic scale Within continuous functional approaches in Coulomb gauge, recent investigations have shown that 
(at least under truncation) there is a direct connection between the Green's functions of Yang-Mills theory and the 
physical string tension that confines quarks [2j. In the Hamiltonian formalism, this relates to both the non-Abelian 
color Coulomb potential |3[, and the temporal Wilson loop [4(. 

On the other hand, the potential that describes the interaction of three quarks is much less studied then the 
potential between a quark-antiquark pair. Lattice simulations have indeed shown that the Yang-Mills interaction 
gives rise to a linearly rising potential between three quarks, but unfortunately a consensus regarding the shape of the 
gluonic strings connecting the three quarks is still missing: either the strings meet at the so-called Fermat point, which 
has minimal distances to the three sources (so-called Y configuration) [5|-l7|, or the qqq potential is simply the sum 
of two-body interactions (so-called A configuration) Q. In the continuum, the only calculations of the three quark 
potential have been performed at perturbative level, within the so-called potential non-relativistic QCD approach 
In Ref. [3, the authors considered the perturbative static potential of three heavy quarks and found that up to 
NLO, the potential is just the sum of the two-body contributions, whereas at NNLO the three-body contributions do 
appear, signaling the importance of three-body interactions for understanding the shape of the string in the infrared 
regime. 

The Faddeev equation [TTj] and its subsequent developments [H, 03 (for an extended review see [I3|) provide a 
general formulation of the relativistic three-body problem. It is a bound state equation (the direct analogue of the 
homogeneous two-body Bethe-Salpeter equation) and it has been efficiently applied in QCD to study baryon states, via 
the Green's functions of the theory. Typically, these studies are performed in Landau gauge and, due to the complexity 
of the equations, they have been mainly restricted to rainbow-ladder truncation, where the kernel is reduced to the 
single exchange of a dressed gluon. Within this approximation and by employing phenomenological Ansatze for 
the Yang-Mills part of the theory, the nucleon and A properties have been studied fla - fla ]. Other simplifications 
include the three-body spectator formalism [l4j, a Salpeter-type equation with instantaneous interaction or the 
diquark correlations [20(. Already from the Bethe-Salpeter studies for mesons, it is known that truncating the kernel 
is not a simple task - the truncation has to be consistent with the symmetry properties of the theory, e.g., the 
axialvector Ward-Takahashi identity must be satisfied. In contemporary studies, the Bethe-Salpeter kernel has been 
considered beyond the rainbow-ladder truncation for meson states, including both vertex corrections |2ll426| as well 
as unquenching effects [22l. [27rl29| and it was found that (apart from meson decay induced by unquenching |28||) the 
rainbow-ladder approximation works surprisingly well. 

For many years, Coulomb gauge studies have been recognized as a promising avenue with which to investigate 
the nonperturbative regime of QCD [301 ] . In this gauge, the Gribov-Zwanziger scenario of confinement becomes 
particularly relevant [3l|-|33|. In this picture, the temporal component of the gluon propagator becomes infrared 
enhanced, providing for a long range confining force, while the transversal spatial component is infrared suppressed, 
thus explaining the absence of the asymptotic states in the spectrum. Coulomb gauge is physical, in the sense that in 
this gauge the system reduces naturally to the physical degrees of freedom (explicitly demonstrated in (33j). Moreover, 
within the first order functional formalism it has been shown that the total charge of the system is conserved and 
vanishing, and the well-known energy divergence problem disappears l34l . Within this approach the Dyson- Schwinger 
equations for the Yang-Mills part of the theory have been derived |35|, [36j , together with the Slavnov- Taylor identities 
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[37j and perturbative results have been provided [38[ . In addition, the quark sector has been also investigated, within 
perturbation theory [39j as well as in the heavy mass limit 2fl. On the lattice, initial calculations for the Yang-Mills 
propagators have also become available [13,0] (see also 45M47| ) . In particular, the results indicate that the temporal 
component of the gluon propagator is largely independent of energy (due to noncovariance, in Coulomb gauge the 
propagators are in general dependent on both the energy and spatial momentum), and it is consistent with a 1/q 4 
behavior in the infrared. Moreover, the spatia l eq ual-time gluon propagator is found to be vanishing in the infrared, 
in agreement with the Gribov's formula (43l . |48| . The lattice calculations sup port the results obtained from the 
variational method to the Hamiltonian approach in Yang-Mills theory [3j, I49l45l| . 

In this paper, we extend a previous investigation of the heavy quark system in Coulomb gauge Q. There, the 
Bethe-Salpeter equation for qq bound states was studied with a heavy quark mass expansion (which underlines the 
Heavy Quark Effective Theory [HQET] [4JJ-|42j) at leading order, and a direct connection between the temporal gluon 
propagator and the string tension was found. Following the same approach, we consider in this work the Faddeev 
equation for three-quark systems in Coulomb gauge, in the symmetric case (i.e., equal spatial separations between 
quarks) and with the inclusion of only two-body interactions, at leading order in the mass expansion. We will use the 
results inspired by the lattice for the Green's functions of the Yang- Mills sector and in addition, we will employ our 
previous findings, in particular that nonperturbatively the temporal quark-gluon vertex remains bare under truncation 
and the kernel of the Bethe-Salpeter equation reduces to the ladder approximation. In this truncated system, we will 
provide an exact solution to the Faddeev equation, from which the confining potential between three quarks emerges, 
and we will show that qqq bound states can only exist for N c — 3 colors, i.e. color singlet baryons. (In [2j it was 
shown that only color-singlet meson and SU(2) qq "baryon" states have finite energy.) 

The organization of this paper is as follows. In Sec. II we briefly review relevant results for the heavy quark systems. 
Starting with the generating functional of Coulomb gauge QCD at leading order in the mass expansion, we review 
the main steps in the derivation of the heavy quark propagator and the temporal quark-gluon vertex needed in this 
work. In Sec. Ill the Faddeev equation for three-quark states is considered. In addition to solving the equation, the 
pole structure of the quark-baryon vertex is analyzed and an interesting similarity with the qq system is discussed. 
Moreover, a direct connection between the temporal gluon propagator and the physical string tension is found. In 
Sec. IV a short summary and the concluding remarks will be presented. Some technical details are given in the 
Appendix. 



II. HEAVY QUARK MASS EXPANSION 

In this section, let us briefly review some relevant results from [2j. The notations and conventions used in this work 
are those established in Refs. pi I35U371 139| . We work in Minkowski space, with the metric g^ — diag(l, — 1). Roman 
letters (i, j . . ) refer to spatial indices and superscripts (a, &,...) stand for color indices in the adjoint representation 
of the gauge group. Unless otherwise specified, the Dirac spinor, flavor and (fundamental) color indices are denoted 
with a common index (a, /?...). Configuration space coordinates may be denoted with subscript (x, y,. . . ) when no 
confusion arises. The Dirac 7-matrices satisfy the Clifford algebra {7^,7"} = 2g ,tv . The notation 7* refers to the 
spatial component, where the minus sign arising from the metric has been explicitly extracted when appropriate. 

The explicit quark contribution to the full QCD generating functional within our conventions, can be written (39| 

Z[x,x] = Jv^cxp^i J d A xq a {x) tr/°D + ij-D - m ^ ^ (a?) | 

x exp \ 1 I d 4 x [x a (x)q a (x) + q a (x)xa(x)) + iSym >■ (2.1) 



In the above, 2?$ generically denotes the functional integration measure over all fields and Sym is the Yang-Mills 
contribution to the generating functional. q a denotes the full quark field, q a is the conjugate (or antiquark) field, 
and Xa, Xa are the corresponding sources. The temporal and spatial components of the covariant derivative (in the 
fundamental color representation) are given by 

D = d -igT a a a (x), 

D = V + tgT a A a (x), (2.2) 

where A and a refer to the spatial and temporal components of the gluon field, respectively. f abc are the structure 
constants of the SU(N C ) group, with the Hermitian generators T a , satisfying [T a ,T b ] = l j abc T c and normalized via 
Tr(T a T b ) = 5 ab /2. For later use we introduce the Casimir factor associated with the quark gap equation: 

N 2 - 1 

C F = — • (2.3) 
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In the following, we briefly sketch the derivation of the quark propagator, in the heavy mass limit. For an extended 
discussion, the reader is referred to the original work [2j. We start by performing the following decomposition of the 
quark field 

q a (x) = e- mxo [h{x) + H(x)] a , h a (x) = e imx ° [P+q(x)] a , H a (x) = e mx ° [P-q{x)] a (2.4) 
(similarly for the antiquark field) and introduce the two components h and H, with the help of the spinor projectors 

P± = i(l±7°). (2.5) 

This corresponds to a particular case of the heavy quark transform underlying HQET with the velocity vector v^ L = 
(1, 0, 0, 0) [40J, but in the functional approach here it can be regarded simply as an arbitrary decomposition that will 
turn out to be very useful in Coulomb gauge (precisely, this will lead to the suppression of the spatial gluon propagator 
in the mass expansion, see below). 

We now insert the decomposition, Eq. (12.41) . into the generating functional Eq. (|2.ip . integrate out the TJ-fields, and 
make an expansion in the heavy quark mass (throughout this work, we will use the standard term "mass expansion" , 
instead of "expansion in the inverse mass"). At leading order, we get the following expression: 



Z[X, x] = »exp M d A xh a (x) [td 0x + gT a a a (x)] a0 hp(x) 



x exp 



d A x [e- mx °x a (x)h a (x) + e* mx °h a (x) Xa (x)] + iS YM } + O (1/m) , (2.6) 

where the temporal component of the covariant derivative D has been written explicitly. In the above, we have 
kept the full quark and antiquark sources (rather than the ones corresponding to the components of the quark field, 
introduced in HQET). This means that we can use the full gap, Bcthc-Salpcter and Faddccv equations of QCD but 
replace the kernels, propagators and vertices and restrict to the leading order in the mass expansion. 

In full QCD (i.e., Coulomb gauge within second order formalism, without the mass expansion and derived from the 
first order formalism results of Ref. [39(), the quark gap equation is given by [dui = d 4 ui / (2tt) 4 ]: 

iW*0 = rf^fc) + / dw {r<%(fc, - U)U - *)wWw)iW( W) -k, k - «,)w£(k - w) 

+ r ?iV( fc ' 00 ~ k ) W mvM T \ q A lSj {^ -k, k - u)Wf Aij {k - w)} (2.7) 

(Waa is the spatial gluon propagator, which will not be regarded here). The gap equation is supplemented by 
the Slavnov- Taylor identity, which follows from the invariance of the action under a Gauss-BRST transform Q ■ In 
Coulomb gauge, this identity reads: 

klv^ikxMM) = i^ qq A*pi(ki,k2,k 3 )r*J(-k 3 ) 



+^qq a s{ki) f|.- C(J (fci + q , fc 3 - q a ; k 2 ) + igT d 
+ \K.ccq( k 2 + qo, h ~ 9o; h) - igT d 
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where ki + k 2 + k3 = 0, go is an arbitrary energy injection scale (arising from the noncovariance of Coulomb gauge [37|), 
Tec is the ghost proper two-point function, Tq.^ cq and T q .^ C q are ghost-quark kernels associated with the Gauss-BRST 
transform. 

Now, as a consequence of the Coulomb gauge decomposition, Eq. (12. 4|) , the part of the generating functional given by 



Eq. (|2.6p corresponding to the tree-level spatial quark gluon vertex T^ q \ is contained within the 0(l/m) contribution 
which is here neglected. Under the further assumption that the pure Yang-Mills vertices may be neglected, the Dyson- 
Schwinger equation for the nonperturbative spatial quark-gluon vertex then furnishes the result that ~ 0(l/m) 
(see for a complete discussion and justification of this truncation). Similarly, the ghost-quark kernels can be 
neglected. Thus, under our truncation scheme, the Slavnov- Taylor identity reduces to 

k° 3 r± qaa0 (fci, fc 2 , k 3 ) = r 5flatf (fci) [igT d ] sp - [zgT d ] aS T qqSp (-k 2 ) + O (1/m) . (2.9) 
This is then inserted into Eq. (|2.7[) , together with the tree- level quark proper two-point function 

r £U fc ) = «W [*o -m} + (1/m) (2.10) 
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and the tree level quark gluon vertex 

4%( fc i' fc 2,fc 3 ) = [gT a ] a0 +0(l/m) (2.11) 

that follow from the generating functional Eq. (|2.6[) . The general form of the nonperturbative temporal gluon propa- 
gator is given by (36J: 

W? a {k) = 5 ab ^-D aa {P). (2.12) 
k 

Lattice results, and also more formal consideration in continuum show that the dressing function D aa has some part 
that is independent of energy 52 1 and moreover, D aa is infrared divergent and likely to behave as 1/k 2 for vanishing 



k 2 (the explicit form of D aa will only be needed in the last step of the calculation). Putting all this together, we find 
the following solution to Eq. (|2.7[) for the heavy quark propagator: 

^ W = jV^l +C)(lH ' (2 - 13) 
with the (implicitly regularized, denoted by "r") constant [du = d 3 uj/(2ir) 3 ] 

2 J r u z 

When solving Eq. (|2.7p . the ordering of the integration is set such that the temporal integral is performed first, 
under the condition that the spatial integral is regularized and finite. Inserting the solution Eq. (|2.13l) into the 
Slavnov- Taylor identity, we find that the temporal quark-gluon vertex remains nonperturbatively bare: 

If^fci.fei, fcs) = [gT a ] a p + O (1/m) . (2.15) 

The propagator Eq. (|2.13|) has a couple of striking features, which have been emphasized in 2]. Firstly, due 
to the mass expansion, we only have a single pole in the complex fco-plane, as opposed to the conventional quark 
propagator, which possesses a pair of simple poles. Hence, it is necessary to explicitly define the Feynman prescription. 
From Eq. (|2.13l) it then follows that the closed quark loops (virtual quark-antiquark pairs) vanish due to the energy 
integration, which implies that the theory is quenched in the heavy mass limit: 

=0. (2.16) 



[fco — m — T r + is] [kg + po — m — I r + is] 



Secondly, the propagator Eq. (|2.13|) is diagonal in the outer product of the fundamental color, flavor and spinor spaces 
- physically this corresponds to the decoupling of the spin from the heavy quark system. In fact, is identical to 

the heavy quark tree-level propagator up to the appearance of the mass term, and this is due to the fact that in 
HQET one uses the sources for the large /i-fields directly, while we retain the sources of the full quark fields. Finally, 
let us emphasize that the position of the pole has no physical meaning since the quark can never be on-shell. The 
poles in the quark propagator are situated at infinity (thanks to X r as the regularization is removed) meaning that 
either one requires infinite energy to create a quark from the vacuum or, if a a hadronic system is considered, only 
the relative energy is important. Indeed, it was shown some time ago [53[ that the divergence of the absolute energy 
has no physical meaning and only the relative energy (derived from the Bethe-Salpeter equation) must be considered. 
It is precisely the cancellation of these divergent constants that distinguishes between physical and unphysical poles. 

We also note that for the antiquark propagator the opposite Feynman prescription is assigned such that the Bethe- 
Salpeter equation for the quark-antiquark states has a physical interpretation of a bound state equation. There, 
the quark and the antiquark do not create a virtual quark-antiquark pair, but a system composed of two separate 
unphysical particles (in the sense that they are not connected by a primitive vertex). Moreover, the Bethe-Salpeter 
kernel reduces to the ladder truncation Q ■ The reason is the cancellation of the so-called crossed box contributions 
(i.e., nonplanar diagrams that contain any combinations of nontrivial interactions allowed within our truncation 
scheme) due to the temporal integration performed over multiple propagators with the same relative sign for the 
Feynman prescription (similar to Eq. (|2.16l) . but in this case the terms originate from internal quark or antiquark 
propagators) . 
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FIG. 1: Faddeev equation for three quark bound states. Solid lines represent the quark propagator, the box represents the 
diquark kernel K and the ellipse represents the Faddeev vertex function with the bound state leg depicted by a triple-line. See 
text for details. 

III. FADDEEV EQUATIONS FOR THREE-QUARK STATES 

Let us now consider the Faddeev equation for three-quark bound states. In this work, we employ only the permuted 
two quark kernels K (which coincide with kernel appearing in the Bethe-Salpeter equation for diquark states) and 
neglect the three-quark irreducible diagrams, i.e., genuine three-body forces. This approximation is also motivated 
by the fact that in the quark-diquark model the binding energy is assumed to be mainly provided by the two-quark 
correlations J54|. In this truncation, the Faddeev equation reads (see also Fig. [1]): 

rVy(pi,P2,P3;-P) = -/ Ak {^.^/(/OW^/^ 

+ K ^'j'(k)W qq0 ,^( P2 + k)W qWl „( P3 - k)r ap „ v ,( Pl , P2 +k,p 3 - k-P) 

+ K aT Y a , [k)W 9riW . (pa + k)W- qqa ' a " (Pi " fc)IV 7 " (Pi -k, P2lP3 + k- P)) (3.1) 

where Pi,P2iP3 are the momenta of the quarks, P — p\ + P2 + Pz is the pole 4-momentum of the bound baryon 
state and T is the so-called quark-baryon Faddeev vertex for the particular bound state under consideration and 
whose indices denote explicitly only its quark content. Due to the fact that in the heavy mass limit the spin degrees 
of freedom decouple from the system, at leading order in the mass expansion the Faddeev baryon amplitude r Q ^ 7 
becomes a Dirac scalar, similar to the heavy quark propagator Eq. (|2.13p . The explicit momentum dependence of the 
kernels K is abbreviated for notational convenience. As in the homogeneous Bethe-Salpeter equation, the integral 
equation depends only parametrically on the total four momentum P. 

As discussed at the end of the previous Section, the kernel K reduces to the ladder approximation (constructed via 
gluon exchange) and it reads 

K 0a .^ 0/ (k)=T^ aa ,w: b Ak)r b ^ =9*T« a> W? a {k)T^ (3.2) 

with the temporal gluon propagator and the temporal quark-gluon vertex given by Eq. (|2.12|) and Eq. (|2.15|) , respec- 
tively. Similar to the Bethe-Salpeter equation for meson bound states, the energy independence of this propagator 
will turn to be crucial in the derivation of the confining potential. 

Let us now have a closer look at the energy dependence of the equation Eq. (|3.ip . In the meson case, since the 
Bethe-Salpeter kernel was energy independent, it was straightforward to show that the Bethe-Salpeter vertex itself 
did not contain an energy dependent part. This observation was then used to calculate the confining potential from 
the qq Bethe-Salpeter equation, via a simple analytical integration over the relative energy variable. Unfortunately 
this approach cannot be extended to baryon states: despite the instantaneous kernel, a relative energy dependence 
still remains and thus one cannot assume an energy-independent Faddeev vertex. Therefore, in order to proceed, we 
make the following separable Ansatz for the Faddeev vertex: 

r Q /3 7 (Pl,P2,P3;-P) = ^a/3 7 ^t(Pl,P2,P3> P)^s(P1,P2,P3; P) (3.3) 

where we have introduced a purely antisymmetric (in the quark legs) color factor ^ (the possible baryon color index 
is omitted) and the symmetric (Dirac scalar) temporal and spatial components T t and r s , respectively. 

Inserting the explicit form of the kernel Eq. (|3.2[) and the quark-baryon vertex Ansatz Eq. (|3.3p , the Faddeev 
equation Eq. (|3.ip can be explicitly written as (for simplicity we drop the label P in the arguments of the vertex 
functions) : 

r a01 (pi,P2,p 3 ) = - g 2 T« T T a p j> TK1 J dkw aa (k)w qq { Pl + k)w qq ( P2 - fc)r t (p? + k ,p° 2 - fc 0) ^)r s (pi + k,p 2 - k,p 3 ) 

+ cyclic permutations, (3-4) 
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where the explicit color structure has been extracted (W^ — S ab W a(7 , Wq qa p = 5 a pWq q ). 
With the help of the Fierz identity for the generators T a 

2 \ T<1 ] a p [T^S-y = SaySgp - -^-6 a p6 S <y, (3.5) 

the color structure can be written as 

K^h'^cl^ = -C B *aPy (3.6) 

with Cb = (N c + l)/2N c , where N c is the number of colors, yet to be identified (i.e., the baryon is not assumed to be 
a color singlet). 

In the next step we perform the Fourier transform for the spatial part of the equation, recalling that the heavy 
quark propagator is only a function of energy. We define the coordinate space vertex function via its Fourier transform 

r.(pl,jS,P3) - JdS 1 d22d2 3 e- l ** 1 ^- g '-^- g *r s @ 1 ,£ 2 ,3 3 ) (3.7) 
(similarly for VF CTCT , as in )2|) such that 

J akW aa {k)T a [fr + = J dx x dx 2 a2 3 e^* 1 ^**^'**W aa (22-2x)T.{a u 3 2 ,2s). (3.8) 

Clearly, the component T s trivially simplifies (as before, we have separated the temporal and spatial integrals, under 
the assumption the spatial integral is regularized and finite) and the equation Eq. (|3.4[l reduces to [^fco = dko/{2w)]: 

T t (p ( lp Q 2 ,p a 3 ) = g 2 C B W a!7 (x2 -xi) J dk W qq (p 1 +ko)Wq q (p^-k )r t {p 1 + ko,p a 2 ~ko,pl)+ cyclic permutations. 

(3.9) 

At this point we make a further simplification, motivated by the symmetry of the three-quark system: we restrict to 
a particular geometry, namely to equal quark separations, i.e. |rj = \x 2 — a?i| = |s?3 — x 2 \ = \xi — x 3 \. By inserting 
the explicit form of the quark propagators, Eq. (|2.13p . we have 

r,«, A* = -tcsw„ m J «o M + h _„ r Lt + + «iH-t-' m -2, + g | + cydk 

(3.10) 

Assuming that the vertex T t is symmetric under permutation of quark legs, an Ansatz that satisfies this equation is: 

r.(rf.Arf) = JC aft _, M+ 1 m+2 , ) + g - < 3 -") 

Since the explicit derivation is rather technical, we only give here the solution and defer the details to the Appendix. 

Notice that in the expression Eq. (I3.11[) there are simple poles (in the energy) present. These poles however do 
not occur for finite energies and cannot be physical. As discussed, this is also the case for the quark propagator. 
Intuitively, when a single heavy quark is pulled apart from the system, the qqq state becomes equivalent (i.e., it has 
the same color quantum numbers) to the qq system in the sense that the remaining two quarks form a diquark which 
for N c = 3 would be a color antitriplet configuration, and hence the physical interpretation of the vertex Eq. p. lip can 
be directly related to the heavy quark propagator Eq. (|2.13[) : the presence of the single pole in Eq. (|3. 1 1|) simply means 
that this cannot have the meaning of physical propagation (this would require a covariant double pole). Moreover, the 
divergent constant I r appearing in the absolute energy does not contradict the physics - the only relevant quantity 
is the relative energy of the three quark system. 

With this Ansatz at hand, we return to the formula Eq. (13.10[) . insert the definitions Eq. (|2.12|) and Eq. (|2 . 14[) 
for W alT (x) and T r , and arrive at the following solution for the bound state energy Pq, in the case of equal quark 
separations: 

Po = 3m + ^g 2 f du [C F - 2C B e^ . (3.12) 
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Since the quarks cannot be prepared as isolated states, the only possibilities for the qqq state are either that the 
system is confined (i.e., the bound state energy P increases with the separation), or the system is physically not 
allowed (i.e., the energy Pq is infinite). From the formula Eq. (|3.12j) and knowing that D aa (uj) is infrared enhanced , it 
is clear that in order to have an infrared confining solution (corresponding to a convergent three-momentum integral), 
the condition 

Cb = -f (3.13) 

must be satisfied. This is fulfilled for N c = 3 colors, implying that ^ a pj = £0^7 and that the baryon is a color 
singlet (confined) bound state of three quarks; otherwise, for N c ^ 3 the energy of the the system is infinite for any 
separation \r\. 

Assuming that in the infrared D aa {uj) = X/uP (as indicated by the lattice data [43TJ471 ] and by the variational 
calculations in the continuum Q), where X is some combination of constants, it is straightforward to perform the 
integration on the right hand side of Eq. (|3.12p . with the result that for large separation \ f\: 

P = 3m+l 9 -^\f\. (3.14) 

This mimics the previous findings for qq systems, namely that there exists a direct connection between the string 
tension and the nonperturbative Yang-Mills Green's functions (at least under truncation). In this case, the standard 
term "string tension" refers to the coefficient of the three-body linear confinement term 033 Also, comparing with 
the result of Ref. 0, we find that the string tension corresponding to the qqq system is 3/2 times that of the qq. To 
our knowledge, no direct comparison between the string tensions of the two systems has been made and hence this 
relation would be interesting to investigate on the lattice. The appearance of three times the quark mass stems from 
the presence of the mass term in the heavy quark propagator Eq. (|2.13p which enters the Faddeev equation, and this 
originates from the fact that in the generating functional Eq. (|2.6[) we have retained the full source terms (in contrast 
with the HQET, where one uses sources for the h- fields directly). 



IV. SUMMARY AND CONCLUSIONS 



In this paper, the Faddeev equation, truncated to include only two-body interactions for three-quark states in 
a symmetric configuration, has been considered. At leading order in the heavy quark mass expansion of Coulomb 
gauge QCD and with the truncation to include only the nonperturbative two-point functions of the Yang-Mills sector 
(and neglect all the pure Yang-Mills vertices and higher order functions) the three-quark confining potential has been 
derived and a direct connection between the temporal gluon propagator and the physical string tension has been 
provided. It was found that, as in the case of qq systems, the bound state energy increases linearly with the distance 
for large separations, and that the coefficient was 3/2 times that of the qq system. 

Due to the absence of the three-body interactions and the restriction to a symmetric configuration, no statement 
about the shape of the confining string (A or Y configuration) can be made. From this point of view, a very interesting 
extension of this work would be to consider general separations between quarks and to explicitly include the Yang- 
Mills vertices in the combined system of Dyson-Schwinger and Slavnov- Taylor identities, together with the Faddeev 
equations, and see if one can extract any information about the shape of the string. Also, by the inclusion of the 
Yang-Mills vertices one can investigate the charge screening mechanism, i.e. study how the value of the string tension 
modifies. 
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Appendix: Temporal component of the quark-baryon vertex 

In this appendix we present the explicit derivation of the energy-dependent part of the Faddeev vertex, Eq. (|3.11[) . 
We start with Eq. (|3.10[) and consider the first of the permutations of the energy integral: 

I = ~ [ ft— T T T fr o , T~-^ Tt ^ + fc °'^ - fc o,^). (A.l) 

J \pl + ko - m - l r + ie\ \p2 - k - m - l r + ie\ 

Using 

1 1 r 1 11 

(A.2) 



[z + a + is] [z + b + is] (b — a) \ z + a + is z + b + ie 

and shifting the integration variables, we find that the integral /, Eq. (|A.1[) . depends only on the momentum p® (and 
implicitly on the bound state energy of the system Pa). Explicitly, it reads (using the symmetry of T t ): 

I= -JH o \ ( , t\i / dk °TT— p cT ^—^^t(P -p° 3 + ko,-k ,p° 3 ). (A.3) 

[Pa - p$ - 2(m + X r )\ J [ka + Po — pi — m — T r + ie\ 

Replacing this in the equation Eq. (I3.10[) . we find: 

T t ( Pl , P2 , P3 ) = -2 9 C B W„mY^ [PQ _ p0 _ 2{m + Ir)] J **> [Po _ p o + kQ _ m _ Ir+l£ y (A.4) 

The form of the equation Eq. (|A.4[) suggests that the function T t can be expressed as a symmetric sum 

TtiplM) = /(P?) + f(P°2) + M), (A.5) 

such that the integral equation for T t (function of three variables) is reduced to an integral equation for the function 
/ (of only one variable). The function f(pf) should be chosen such that the integral on the right hand side of the 
equation Eq. (|A.4| generates a factor proportional to \Pq — p® — 2(m + X r )] , to cancel the corresponding factor in the 
denominator. To examine this possibility, we impose the following condition: 

dkQ f(Po-P° + ka) + f(-k ) + M) = . m 



[P Q - p° - 2(m + l r )} J [Pa-pV + ka-m-lr + ie] ~ P - 3(m + 2 r )" 

where a is a (dimensionless) positive constant which remains to be determined. Rearranging the terms to factorizc 
the function f(ko), the above equation can be rewritten as 



dk f(ko) 



1 



= (2a - l)Pg - 2ap\ + (3 - 4a)(m + l r ) „ 
ka-m~l r + is ' Pa-p Q l -ka-m,-l r +ie\ [ % ' 2[P - 3(m + l r )} J[Pl) ' 

Then the most obvious Ansatz for the function / is 

^ k °^ = (2a - 1)P Q - 2ak Q + (3 - 4a)(m+l r ) + ie 

such that on the right hand side of the equation Eq. (|A.7|) the numerator is cancelled by f(pi). The next step is to 
complete the integration on the on the left hand side, which gives (note that the e prescription is chosen such that 
only the first term in the bracket survives - the integration must not give rise to any new terms containing the energy 

1 i 

dk of( k a)z ttz — = -TT, TT5 Z~R — TT? ( A ' 9 ) 

ka — m — X r + ie (2a — l)Po + (3 — 6a) (m + X r ) 

It is then straightforward to compare Eq. (|A.7[) and Eq. (IA.9I) and find that the equality is satisfied for a = 3/2, 
leading to the expression for the vertex F t used in the text. 
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